Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



O 

o 

(N 



> 
m 
in 

m 
o 
\o 
o 

■ 

Oh; 

6 



MORPHOLOGY OF SPITZER 24/iM-DETECTED GALAXIES IN THE UDF: THE LINKS BETWEEN THE 

STAR FORMATION AND GALAXY MORPHOLOGY 

Y. Shi, G. H. Rieke, C. Papovich, P. G. Perez-Gonzalez, E. Le Floc'h 

Draft version February 5, 2008 

ABSTRACT 

We have studied the morphologies of luminous infrared galaxies (LIRGs; L/ii;(8-1000^m) > IO^^Lq) 
at 0.3 < z < 1.4 in the HST Ultradeep Field (UDF) by calculating concentration and asymmetry 
indices and comparing the results with similar calculations for: (1) galaxies at similar redshift that 
are less infrared-active; and (2) local LIRGs. We find that the high-rcdshift samples arc dominated 
by galaxies with concentrations similar to local late-type disk galaxies; however, they are significantly 
more asymmetric than most local galaxies but are similar in both regards to local LIRGs. On aver- 
age, the high-redshift infrared-active galaxies are slightly more asymmetric than the less-active ones, 
although they do include a significantly higher portion of highly asymmetric (merging?) systems and 
a lower portion of more concentrated, symmetric ones. The morphological similarity of infrared-active 
and typical infrared-inactive galaxies at high-redshift suggests that they may be from the same parent 
population, but are in different stages of an episodic star formation process. The similarity between 
high-redshift and local LIRGs suggests that a certain level of asymmetry is generally associated with 
LIRG-level activity. 

Subject headings: galaxies: interactions — infrared: galaxies 



1. INTRODUCTION 

Two fundamental aspects of galaxy formation and evo- 
lution are mass assembly and morphology transforma- 
tion. While mass assembly is often accompanied by 
strong episodes of star formation, the morphology of a 
galaxy refiects its dynamical history and evolution. For 
example, different Hubble types are associated with dif- 
ferent motion patterns of stars and gas and different his- 
tories of star formation. The interplay between these 
two processes may play a crucial role in galaxy for- 
mation and evolution: the dynamics, indicated by the 
morphology, drives the star formation, while the deple- 
tion of gas from star formation as well as feedback from 
young stars through supernovae and stellar winds may 
re-shape the galaxy morphology. Thus, in the local Uni- 
verse, later-type galaxies are gas-rich a nd the sites of el - 
evated star formation (for a review, see iKennicuttll 199^ . 
while the more extreme star-forming galaxies often 
show peculiar morphologies, resulting from collisions and 
mergers (e.g. Soifcr ct al. 1984; Cutri fc Mc Alary 1985 
lenders et al. 1988; An nus et al.lll987HKleinmann et all 
ri98?^.BaTton et a,1...2nn(f ^ 

The comoving star formation densit y rises by an or- 
der of magnitude from z = t o z ^ 1 (Lillv et al. '1995J 
lEUis ct al. 1996; Stcide l et al] IT999: Hopkins 2004). In 
the high-redshift universe, the high rate of star formation 
seems to be associated with a large fraction of peculiar 
and morphologically disturbed galaxies. Deep surveys, 
especially those carried out by the Hubble Space Tele- 
scope (HST), reveal that the fraction of irregular and 
interacting galaxies incre ases dramatically to ward high 
redshift, from around 4% l|Marzke et al.l 19981) in the lo- 
cal u niverse to more than 80% at z >2 ( Conselice et al.l 
l2005j) . Even the high-redshift galaxies classified as tra- 
ditional Hubble types still show many peculiar features. 
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with tid al or disturbed structures and distorted disks. 
However. iMelbourne et al. (2005) study the optical mor- 
phology evolution of LIRGs, finding that peculiar-to- 
spiral LIRG fraction increases at lower redshift. Our 
goal in this paper is to probe how the star formation 
and galaxy morphology arc linked out to z^l. 
The .Hu bble (1926, 19 3 6) galaxy classification a s mod - 



Sandage fc Tammannl lll987tl . 



van den Berghl l)1960albD and 



ified bv iSandaed Il96 
Sanda.ge B'^dO 
ie Vaucou leural l)1959r ). has achieved remarkable success 
in categorizing local galaxies with normal morphology 
(for a review, see ^ Roberts & Hayncs 1994). However, 
the Hubble scheme has a single class for peculiar mor- 
phologies and thus is unable to reflect the vast varia- 
tions within this category. It is therefore inadequate to 
classify the morphologies of high-redshift galaxies using 
the local Hubble sequence, because a significant frac- 
tion of the population are irregular gala xies or nor- 
mal galaxies but with perturbed features ([Driy ^r ^1^ ^1 



1995. 11993 lAbraham et aIjri99fialH: IBrinchma,nn et al 
19981 iCassata et al.l l2005j) . A different classification 



method, called the CA system, is particularly use- 
ful for peculiar galaxies at high redshift. It is based 
on morph ological parameters: the light concentration 
index. C JOkamura et all 119841: lAbraham et al] Il994t 
iBershadv et alJl200rt 'Conselice"2003V and the asvmme- 
try index. A, first introduced by Abraham ct al. (199^) 
and developed bv IConsehce et al.l l|2000j) . Its primary 
virtues are that it is objective and can be applied consis- 
tently to low resolution and low signal-to-noise images. 
The success in distinguishing between different Hubble 
types; the robust measurements of parameters over a 
large range of redshift range of galaxy types especially 
quantifying the irregularity of peculiar galaxies; and the 
strong relations between the parameters and physical 
properties such as color, all make the CA system use- 
ful for classifying galaxies found in deep fields. 
Out to z ~ 1.5, the star formation density is dominated 
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by LIRGs l|Le Floc'h et alJ 120051: iPerez-Gonzalez et alJ 
HqOS) . Studies of the morphologies of these high-redshift 
LIRGs are necessary to understand how morphologies 
and star formation are coupled to drive galaxy evolution. 
In this pa per, we ut ilize a Multiband Imaging Photome- 
ter fMIPS:l Riekeet al., 2004.^ deep survey of the HST Ul- 
tradeep Field (UDF). The extremely deep HST images 
in the UDF allow us to apply the CA system to analyze 
the optical morphologies of the LIRGs in the MIPS sur- 
vey even including faint features with low optical surface 
brightness. 

The morphologies of galaxies in the UDF have 
been recently studied by lElmegreen et alJ |)200,'tH who 
mainly focus on the dis tribution of galaxy types and 
iMenanteau et all 1)2005(1 who study the contribution of 
various types of galaxies to the cosmic star formation 
rate (SFR). The detailed analysis of the morphologies of 
the high-redshift LIRGs and the comparisons to other 
MlPS-non-detected galaxies and local LIRGs in this pa- 
per complement these studies and enable us to gain new 
insight to the following questions: How do the intensely 
star-forming galaxies at z ~ 1 relate to other galaxies 
at the same epoch? How are star formation and galaxy 
morphologies coupled to drive galaxy evolution? 

In this paper, we measure the concentration and asym- 
metry indices from the HST images to quantify the mor- 
phologies of the sample of high-redshift LIRGs in the 
UDF, of a comparision sample of MlPS-non-detected 
galaxies in the UDF, and of a local sample of LIRGs. In 
§ 2, we summarize the observations, data reduction, and 
band merging. In Section § 3, we describe the selection 
of the samples and the measurement of the structural 
parameters. In § 4, we present the morphological charac- 
teristics of the galaxies. In Section § 5, we discuss the im- 
plications for high-redshift galaxy evolution based on the 
results from the preceding section. In § 6, we present our 
conclusions. Throughout this paper, we assume Ho=70 
km s"^ Mpc"\ r2o=0.3 and fiA=0.7. We use AB magni- 
tudes throughout, where mAB = 23.9 — 2.5 \og{f^/l fiJy). 
We denote magnitudes from the HST ACS passbands 
F435W, F606W, F775W, and F850LP as Vgoe, «775, 
and 2850, respectively. 

2. THE DATA 

2.1. Spitzer Data 

Our 24 fiiji MIPS l|Rieke et alJl2004|) observations of 
the UDF are part of a scan map covering a larger 
area centered on the Chandra Deep Field South (CDFS; 
iGiacconi et all 1200^ . The images were processed with 
the MIPS instrument team Data Analysis Tool (DAT; 
[Gordon et al. 2005). Descriptions of 24 ^m source detec - 
tion and photometry are given bv lPapovich et alJl|200^ . 
Table 1 lists the sources detected at 24^m with unam- 
biguous optical counterparts that are the primary sub- 
ject of this paper. The CDFS was also observed un- 
der GT O time with the Infrared Array Camera (IRAC) 
(|Fazio et al.ll2004() at 3.6, 4.5, 5.8 and 8.0 ^iul. Source 
detection and pho tometry were done with SExtractor 
lIBertin fc Arnoutlfl996') following similar procedures as 
in lHuang et alJ l|2004^) . 

2.2. HST Data 

The UDF is a public HST survey to image a single 
Advanced Camera for Surveys (ACS) wide field cam- 



era (WFC) field (11.5 arcmin^) in 4 broad-band fil- 
ters ( F435W: F606W; F775W; F850LP) l|Beckwith et alJ 
I2005j) . The WFC has a resolution of ~0.12", correspond- 
ing to 540 pc at z=0.3 and 1 kpc at z=1.0. For our 
analysis we use the reduced UDF data vl.O made pub- 
lic by the Space Telescope Science Institute (STScI) on 
09 March 2004. The z-band-based catalog is used for 
optical identification of MIPS sources. A total of 7016 
objects detected at z-band gives a source density of 0.169 
arcsec"^. The 10-cr limiting magnitudes at B-band, V- 
band, z-band and z-band are 28.7, 29.0, 29.0 and 28.4, 
respectively, in an aperture of 0.2 arcsec^, correspond- 
ing to 10-cr surface brightness of 27.82, 28.12, 28.12 and 
27.52 arcsec"^, respectively. 

We also use Near-Infrared Camera and Multi-Object 
Spectrometer (NICMOS) UDF Treasury Observations 
for the photometry data for our photo-z code. NIC- 
MOS UDF Observations use camera 3 with a resolution 
of 0.2"/pixel (1.6 kpc at z=l). Due to the small field 
of NICMOS camera, the NICMOS UDF only covers a 
subsection (5.76 arcmin^) of the optica l UDF. The data 
reduct ion and photometry are given bv lThompson et alJ 
((2005|) . The catalog contains 1293 objects, giving a 
source density of 0.0624 arcsec"^. The 5-a limiting AB 
magnitude is 27.7 at 1.1 and 1.6 /im in a 0.6" diameter 
aperture, corresponding to a surface brightness of 27.01 
arcsec"^. 

2.3. Results and Band Merging 

We identified 52 sources with 24 /im flux density 
greater than 0.06 mJy ( a t which level the com plete- 
ness is ^ 50% in the UDF lIPaDovich et al.ll200l ). We 
searched for ACS counterparts on the z-band image of 
the UDF within a 2.5"-radius circular region surround- 
ing each 24/im source. The choice of this search radius 
is motivated by the large FWHM of the MIPS 24 /im 
PSF (~6") and the positional accuracy at 24 fim ( ^ 
0.5" rms). 2.5" (10 kpc at z=0.3 and 20 kpc at z=1.0) 
is large enough to account for the possible physical shift 
between the locations of infrared emission and the opti- 
cal emission observed in local interacting galaxies (e.g. 
ILe Floc'h efsIMM) . 

Multiple optical sources are present in the search cir- 
cle for 41 of 52 MIPS sources. The IRAC images were 
used to identify the most probable optical counterpart for 
these 41 MIPS sources. In the UDF, there are a total of 
312 objects in the IRAC 4.5 /im catalog with flux density 
greater than 1.6 /xJy at 50% completeness. 23 of the 41 
MIPS sources have such an IRAC counterpart within a 
1.0" radius. The probability of random spatial superpo- 
sition of IRAC sources in a 1.0"-radius circular region is 
2.3 %, which indicates a high probability that the IRAC 
counterparts are physically related to the MIPS sources. 
We also visually inspected the IRAC image of each source 
to check if there is nearby IRAC source which might 
contaminate the 24 /im emission. MIPS-15950, which 
is contaminated in this way, is excluded. The search ra- 
dius to match optical counterparts for IRAC sources is 
1.5", motivated by similar reasons to those for the MIPS 
sources. 15 of 22 IRAC sources have only one optical 
source within the search radius, which should be the op- 
tical counterpart of the corresponding MIPS source. For 
the remaining 25 MIPS sources with possible multiple 
optical counterparts, if an optical object is brighter than 
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other optical objects within the same search circle by 
4 magnitudes, this optical object was taken to be the 
optical counterpart of the MIPS source. This associ- 
ation is justified because IR-luminous gala xies tend to 
be associated with optically bright galaxies l|Mann et al.l 
IT997.,') . On the other hand, the probability of random 
spatial superposition is an order of magnitude larger for 
4-magnitude fainter objects. Based on the above series 
of steps for deconfusion of multiple objects in the 24 /xm 
beam, 16 of 52 MIPS sources may still have multiple op- 
tical counterparts. We do not use these 16 sources in our 
analysis to avoid their affecting the distribution of the 
morphologies of the MlPS-detected sources. Our final 
sample therefore consists of 35 MlPS-detected galaxies 
that have unambiguous identifications with ACS images. 

An alternative method for evaluating interactions, 
"pairs statistics", would be strongly affected by the 
handling of sources with multiple possible counterparts. 
However, the CA analysis employed here measures the 
structure of the central image and hence to first order is 
independent of such effects. As a test, we recomputed the 
classification parameters including the objects with mul- 
tiple identifications and using the brightest object as the 
optical counterpart. We found that the basic results were 
unchanged; the mean asymmetry of the MlPS-detected 
objects got slightly smaller but by only 0.03. 

Eighteen of our sample of 35 MIPS sources have 
spectroscopic red shifts based on the y LT/F0RS2 spec- 
troscopy catalog l|Vanzella et alJl2Q0,'iD and VIMOS VLT 
Deep Survey (VVDS) vl.O catalog (! Le Fevre et all200l . 
For the remaining 17 sources, we use photometric red- 
shifts. The publicly available COMBO-17 survey of the 
CDFS provides p hotometric redsh ifts by means of 17- 
band photometry ( W olf et aril2004H . The typical photo- 
metric redshift accuracy is 6z/{l + z) < 0.02 for galaxies 
with mji < 22 and reaches 0.1 at toa, ~24, allowing the 
rest- frame absolute magnitude to be accurate within 0.1 
for mji < 22 and 0.5 mag for mn ~24. We searched for 
COMBO-17 counterparts for each ACS source within a 
radius of 0.45 arcsec. Given a total of 63,501 COMBO- 
17 objects in a field of view 31.5'x 30', the probabil- 
ity of random spatial superposition in a 0.45"-radius 
circular region is only 1.1%. We visually inspected 
the ACS image of each source to assure we matched 
the right COMBO-17 counterpart. Seven MIPS sources 
have X- ray emi ssion based on the X-ray 1 MS catalog 
(j^iacconi et ^ '2002^ and photometric redshifts given 
bv W. Zheng et al. ( 2004 ). For the ACS counterpart of 
any MIPS source without redshif ts from these sources, 
we used our own photo-z code fsee lPerez-Gonzalez et al.l 
EaOS) based on ACS, NICMOS and IRAC observations. 
We also used our code to compute the redshift for any 
ACS source with a COMBO-17 counterpart with niR > 
22. The main advantage of this photo-z code is it uses 
photometry over a longer wavelength baseline (includ- 
ing the 1.6 Ilia bump) to obtain redshifts for galaxies 
up to Z'^3, with accuracy of Sz/{1 -I- z) < 0.1. We in- 
spected the fitting individually to assure our code gave a 
robust measurement of redshifts. Except for MIPS-5106, 
sources with multiple redshifts have consistent measure- 
ments of redshift. For this source, we use our redshift 
measurement since the 1.6 fim stellar bump is obviously 
detected at the IRAC bands in the SED. 

The K-corrections for the optical four-band data are 



computed using codes based on the template spectra 
of galaxies in the Sloan Digital Sky Survey (SDSS; 
iBlanton et al]l200l . 

We used the ACS image that is closest to rest-frame 
B band to compute the structural parameters for the 
galaxies. Therefore the maximum redshift in our study 
is 1.4, beyond which the observed z-band is shorter than 
the rest-frame B-band. The advantages of classifica- 
tion using rest-frame optical bands and the importance 
of morphological K-cor recti ons are discus sed further in 
iWindhorst et all ()2002D and lPapovich et al.. (^003). 

To obtain accurate total far-infrared luminosities, it 
is necessary to have photometry over the whole range 
of infrared wavelengths. In this study, however it is 
enough to have a coarse estimate of the total infrared 
luminosity. We used the template spectra of starburst 
galaxies (with the current SFR much larger than the 
past-averaged SFR) pubhshed bv iLagache et al.l (|2004fl 
to get the 8-1000 /xm total infrared luminosity from the 
observed 24 fim flux density. Due to the high sensitivity 
of MIPS, some MlPS-detected galaxies at low redshift are 
normal galaxies, not starbursts. For them, we used the 
appropriate template for galaxies with total IR emission 
fainter than 10^° Lq (where the threshold is computed 
from the starburst template). The average conversion 
factor from i'L^{2A fim) to total infrared luminosity is 

20 with a dispersion of up to 0.2 dex, indicating that 
our conver sion should be accurate typically within a fac- 
tor 2 or 3 (|Par)ovich fc Bel]ll200llLe Floc'h et alJl2b051 
which is enough for the analysis in this paper. 

3. METHODOLOGY 

3.1. Sample Definition 

We only computed the structural parameters for 34 of 
the 35 MIPS sources that have signal-to- noise (S/N) ra- 
tios better than 100 and half light radius r^o > 5 pixels 
at z band. High signal-to-noise ratio (SNR) and large 
size are req uired for robust measures of the asyrnm etry 
parameters (|Consehce et al.ll2000l iBershadv etalE oOOl . 
The deep exposure in the UDF thus allowed us to com- 
pute reliably structural parameters for a high fraction 
(97%) of MlPS-detected galaxies. 

The members of our control sample are defined to be 
undetected at 24 fim and (1) to lie 5 " away from all 
MIPS sources; (2) to have half light radius rso > 5 pix- 
els at z band; (3) to have SNR better than 100 at z 
band; and (4) to have a COMBO-17 photometric redshift 
smaller than 1.4. The last condit ion is necessary for K- 
corrections and requires niz < 25 (jWolf et alJl2004|) . The 
final sample is composed of 252 normal galaxies, repre- 
senting 80% of the galaxies undetected at 24/im meeting 
only the first three requirements. 

Our main goal is to investigate the morphologies of 
high-redshift LIRGs. We identified a galaxy with Ljjf 
> IO^^Lq as a LIRG; the final sample of LIRGs has 

21 sources. The deep exposure of the UDF allows us 
to detect their optical counterparts down to low surface 
brightness, removing potential biases from less-deep clas- 
sification imagery. The 10-ct surface brightness at z-band 
is 27.52 arcsec"^, which corresponds to a rest-frame B- 
band surface brightness of 24.5 for a fiducial galaxy at 
z=l. This corresponds well to the typical surface bright- 
ness level reached in morphological studies of local galax- 
ies (although at lower physical resolution). It is also close 
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to the central surface brightness of a local low-surface- 
brightness (LSB) galaxy. In addition, the LIRG sample 
is sufficiently large for reliable statistics, as can be seen 
from the results in Table 2. 

We also created a sample of 5 MlPS-detected but lower 
luminosity non-LIRGs and a control sample of 137 MIPS- 
non-detected galaxies. The sample of MlPS-detected 
non-LIRGs is composed of galaxies with Lm < 10^^ Lq 
and Mb < -18.5. The sample of MlPS-non-detected 
galaxies is derived from the sample of 252 normal galax- 
ies with the additional constraints that Mb < —18.5 and 
that the redshift be greater than the minimum redshift 
of the LIRGs. The cut-off magnitude Mb = —18.5 is ap- 
proximately the rest-frame B magnitude of the faintest 
LIRG. 

Figure shows the redshift distributions of the three 
samples and the infrared luminosity at 24 /im flux den- 
sity of 0.06 mJy as a function of redshift. Galaxies of 
all three samples mainly fall at z>0.3. The sample of 
LIRGs and MlPS-non-detected galaxies have almost the 
same redshift distributions while non-LIRGs mainly fall 
at lower redshift due to the detection limits. In the fol- 
lowing analysis, we mainly discuss the morphologies of 
the LIRGs and comparisons to the MlPS-non-detected 
galaxies, since the MlPS-detected non-LIRGs have a dif- 
ferent redshift range and are few in number. 

3.2. Quantitative Morphology Classification: 
Concentration and Asymmetry 

Before measuring the concentration and asymmetry 
value for each galaxy, some image processing was re- 
quired. We identified and masked contaminating sources 
in the region of the galaxy using 2;-band segmentation 
maps provided with the UDF data release. We replaced 
the mask region with noise to compute the morphology 
parameters. 

We measured the con centration index usin g the 
methodology described in iBershadv et al.l 1)2000(1 . To 

study the morphologies in a fixed-size aperture at the 
physical scale of each galaxy, the dimensionless param- 
eter ?7 as defined by ri{r) = I{r)/ < I{r) > (Petrosian 
11976(1 is used to define the total size of the galaxy, where 
/(r) is the surface brightness at radius r and < I{r) > is 
the mean surface brightness within radius r. The appar- 
ent total magnitudes are then integrated over an aperture 
equivalent to the total size of the galaxy defined as twice 
the radius r(ry = 0.2). Based on the measured curve of 
growth, the concentration index is defined by 

C = 51og^, (1) 

where rgo and r2o are the radii that enclose 80% and 20% 
of the total light (.Kent.1985) , respectively. 

The full descriptior i of the algorithm to compute the 
asymmetry index is in iConselice et"alll|2000|l . The asym- 
metry parameter is defined as: 

A^mm{ =rr7-^ ) - mm( =r— ), (2) 

where the sums are over the total size of the galaxy as 
defined above. In the above equation, the first term rep- 
resents the asymmetry of the galaxy and the second one 
corrects the effect of noise on the asymmetry measure- 
ment, where Iq and /igo are the intensity of each pixel 



of the image and of the image rotated by 180°, respec- 
tively, and Bq and Siso are the intensity of each pixel of 
the background region and its 180°-rotated counterpart, 
respectively. We extracted a 100x100 pixel background 
region without any contamination and scaled the asym- 
metry of this region by the relative size of the galaxy to 
this background region area. The asymmetry minimum 
was obtained by searching the minimum of asymmetries 
at rotation centers a round the initial galaxy center. See 
IConselice et al.l ((2000) for the detailed computation algo- 
rithm to locate the rotation center producing the asym- 
metry minimum. 

The uncertainty in the concentration measurem ent is 
mainly due to the spatial resolution. Bcrshadv et alJ 
((2000fl show that the uncertainty of this parameter is 
around 0.2 for >5 pixel. For asymmetry mea- 
surements, the noise of the image is the main effect. 
IConselice et al.l 1(2000(1 show that asymmetry measure- 
ments have errors around 0.02 for S/N>500 and around 
0.05 for 100<S/N<300. We tested our programs by com- 
puting the concentration and a symmetry parameters of 
20 galaxies in iConselic"^ l(2003ll with asymmetry span- 
ning from 0.01 to 0.4 and concentration from 2.5 to 4.5. 
The average differences in concentration and asymmetry 
parameters are 0.2 and 0.02, respectively. For galaxies 
with SNR better than 100 in our sample, we estimate 
the typical errors of the concentration and asymmetry 
indices are around 0.3 and 0.04, respectively, which are 
adequate for our study. 

4. RESULTS 

4.1. Morphologies of a Representative Sample of LIRCs 
at 0.3<z<14 

An important aspect of this study is that we can deter- 
mine morphological properties for a representative sam- 
ple of high-z LIRGs. In the UDF, 23 of 34 MlPS-detected 
galaxies at 0.3<z<1.4 with Mb < -18.5, S/N > 100 and 
> 5.0 are LIRGs. Compared with 137 MlPS-non- 
detected galaxies with the same constraints as above, 
we find that ~15% of the galaxies with Mb < —18.5, 
S/N > 100 and r^o > 5.0 are LIRGs at z~l. This frac- 
tion may be underestimated due to the non-detection 
of LIRGs near Lju ~ 10^^ Lq and also may be over- 
estimated as a result of missing non-LIRGs near Mb = 
— 18.5 at high redshift. However, the portion of LIRGs in 
our sample is similar to those found by others at similar 
redshift. For example, the fraction of LIRGs at z>0.4 
in the Canada France Redshift Survey (CFRS) is about 
16% without correction for the non-detection of galax- 
ies near the cut-off limits fHammer et al.ll200 5f). An- 
other test is to use the in frared and optical luminosity 
functions. iLe Floc'h et all l(2005tl recently obtained the 
infrared luminosity function at 0<z<1.2 based on the 
MlPS/Spitzer 24 /xm sources located in the CDFS. Us- 
ing L*jjj = 1.9 X IO^^Lq, aiR = 1.23 and ajR = 0.72 at 
1.0<z<1.2, we estimate 70% of the LIRGs at 0.9<z<1.4 
are below 0.06 mJy. Using the rest-frame i?-band lu- 
mino sity function of the COMBO- 17 survey l(Wolf et alJ 
12003(1 . we estimate 65% of the galaxies with Mb < —18.5 
are below tub = 25. Correcting the fraction of non- 
detected galaxies at high redshift, we find that the frac- 
tion of LIRGs is almost the same (~18%) as in our study. 
At z< 0.9 where both LIRGs and Mb < —18.5 galaxies 
are detected, we obtain a fraction ~15% of LIRGs for a 
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total of 67 galaxies. We conclude from all these argu- 
ments that the fraction of LIRGs for Mb < —18.5 is ~ 
15% at 0.3<z<1.4, the same as in our study. This frac- 
tion is much larger than that in the local uni verse, which 
is 0.5 % for galaxies with itot > 10^" Lq ijSoifer et al.1 

[Mi). 

Table 1 lists the LIRGs detected in the UDF along with 
their concentration and asymmetry indices, C and A. 
MIPS-4641, MIPS-5117 and MIPS-15236 have X-ray de- 
tections (Giacconi ct al. 2002) and are classified as type 
1 AGN bv .W. Z heng ct al. ( 2004). These three sources 
are excluded in our study of morphologies. The other 
two MIPS s ources detected in the X-ray are classified as 
galaxies bv IW. Zheng et aP l|2004|) and are kept in this 
study. 

Table 2 lists the statistical results for LIRGs, non- 
LIRGs and the MlPS-non-detected sample in the inter- 
mediate redshift bin 0.3<z<0.9, in the high redshift bin 
0.9<z<1.4, and for the local (low redshift) LIRGs. At 
0.3<z<0.9, the 24 detection limit (0.06 mJy) and the 
B-band magnitude limit (-18.5) yield complete samples 
of LIRGs and MlPS-non-detected galaxies, respectively. 
At 0.9<z<1.4, both the LIRGs and MlPS-non-detected 
samples are incomplete. We mainly discuss our result 
based on the complete sample. However, the MlPS-non- 
detected sample at 0.9<z<1.4 is still dominated by non- 
LIRGs. Given 9 LIRGs, 88 MlPS-non-detected objects 
and 70% incompleteness in IR-detection at 0.9<z<1.4, 
the fraction of LIRGs in 88 MlPS-non-detected objects 
is <25%. Therefore, the comparision of morphologies 
between LIRGs and MlPS-non-detected galaxies at this 
redshift range is still helpful in our understanding the dif- 
ference in morphologies between infrared-active galaxies 
and infrared-inactive galaxies. 

Figure [3 shows the rest-frame _B-band images for 
LIRGs used for our morphological study. Figure |21 plots 
the galaxies along with our other two high-z samples on 
the C, A plane. The solid curves show the division of this 
plane into three regions mainly populated by early-type, 
intermediate- type and late- type galaxies, based on clas- 
sifications of local galaxies (c.f. Abraham ct al. 1996a 0: 
iConselice et al.M2000() . Figure |2| shows that the high- 
redshift LIRGs are predominantly lat e-type galaxies (as 
found previously bv lBell et alJl2005aj) with asymmetric 
structures, including long tidal tails and distorted disks. 
Some of them have very low surface brightness, such 
as MIPS-5106 with rest-frame B-band surface brightness 
/iBo = 22.4 within the Petrosian radius. 

We define a galaxy with A>0.2 as an asymmetric 
galaxy in this study. This is motivated by the fact that 
most (^^90%) of local normal galaxies including elliptical 
galaxies, early-type (Sa-Sb) and la te-type spiral galaxies 
(Sc-Sd), have asymmetry A<0.2 ijConsehce .2003.) . As 
hsted in Table 2, 90l:3!^% of the LIRGs in the UDF are 
late-type galaxies and 80±10% of the LIRGs are asym- 
metric (A>0.2) objects. 

The classification of the high-z galaxies using the asym- 
metry and concentration indices is generally consistent 
with the conventional visual classifications. Figure El 
shows that only two sources, MIPS-5115 and MIPS- 
13901 are classified visually differently from the results 
using the structural parameters. In this paper, we make 
no use of the conventional classifications (other than to 



place some results in context), basing all of our conclu- 
sions on the CA analysis. Thus, the small number of 
classification discrepancies has no effect on our conclu- 
sions. 

4.2. Comparision of LIRG Morphologies to those of 
MlPS-non- detected Galaxies 

In this Section, we compare the morphologies of LIRGs 
in the UDF to those of MlPS-non-dctectcd galaxies. Ta- 
ble 2 shows that the correlations at high redshift between 
infrared luminosity and morphology are weak. On aver- 
age, high-z LIRGs are more asymmetric than the MlPS- 
non-detected galaxies, ch aracterized by a higher frac- 
tion oi A> 0.35 (merging: |C_onselice 2003) galaxies and 
fewer symmetric objects, defined as A<0.2. However, a 
number of MlPS-non-detected galaxies show asymmet- 
ric structures and/or are currently involved in merging 
activity. In the following we compare the morphological 
distributions in terms of the fraction of late- type galaxies 
and the distribution of asymmetry parameters. 

4.2.1. Distributions of the Asymmetry and 
Concentration Parameters 

Table 2 lists the average asymmetries and the standard 
deviations of the means. The LIRGs have asymmetry 
around 0.34, somewhat higher than A=0.22 of the MlPS- 
non-detected galaxies. The K-S test shows the probabil- 
ities that LIRGs and MlPS-non-detected galaxies have 
the same distribution of asymmetries are 1% and <0.1% 
for the intermediate-redshift and high-redshift bins, re- 
spect ively. Th is seems consistent wi th the suggestion 
(e.g. iLarson fc Tinslcv 1978: Sande rs fc Mira bel 199^ 
iBarton et afWoOQ.) that non-symmetric structures, such 
as tidal distortions, companions, and distorted disks, 
play an important role in enhancing star formation in 
galaxies. Figure 13 indicates that the MlPS-non-detected 
galaxies have a wide distribution of morphologies, from 
classical Hubble types to highly asymmetric, most likely 
merging systems. The fraction of symmetric (A<0.2) 
galaxies in this sample is high (~ 50 %). In comparison, 
the LIRGs virtually lack symmetric objects (~10%). 

As listed in Table 2, the sample of LIRGs shows a 
slightly higher fraction of late-type galaxies than the 
sample of MlPS-non-detected galaxies, although the val- 
ues are consistent within statistical errors. The classifi- 
cation calibration method in this comparison is based o n 
visually classi fied local galaxies |Conselice et al.l 12000'). 
Alternatively, iGassata et al.l l|200 5h measured the rest- 
frame B morphology for the K-band selected high- 
redshift galaxies and used the boundary {A <0.2 and 
C >2.9) to define the early-type galaxies. Based on their 
method, we find that the late-type fraction of LIRGs 
is still a little higher at 0.3<z<1.5, i.e., the 24/xm- 
selected sample is biased a little toward late-type galax- 
ies. However, much of this relatively minor difference 
can be explained by the larger portion of symmetric 
early- and intermediate-type galaxies in the MlPS-non- 
detected sample. The differences in morphologies be- 
tween the late-type galaxies of these two samples are 
therefore minor. 

On the other hand, the galaxies of both samples 
are more asymmetric than the local Hubble-type galax- 
ies and even slightly more asymmetric than the non- 
starburst local dwarf irregular galaxies, which have 
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asymmetry of 0.17±0.10 ijConselicd 1200^ . The K-S 
test shows that the probabihty that LIRGs and local 
irregular galaxies have the same distribution of asym- 
metry is <0.1%, while such a probability for MIPS- non- 
detected g alaxies at 0.3<z< 0.9 and local irregular galax- 
ies is 27%. iConselice et alJ 12000) show that low spatial 
resolution affects the asymmetry measurement by mak- 
ing the galaxy appear more symmetric. The spatial res- 
olutions are 50 pc/pixel for the local sample of Consclicc 
(|2003f) and 300 pc/pixel for our galaxies at z=0.5. This 
reinforces our conclusion that the galaxies in the UDF 
are more asymmetric than the local galaxies. 

Although the two high-z samples show somewhat dif- 
ferent levels of asymmetry, their concentrations are al- 
most the same, characterized by values that are compara- 
ble to the local late-type disk-galaxies (Sc-Sd; Conselice 
|2003). The K-S test shows that the LIRG sample and 
the MlPS-non-detected sample have the same distribu- 
tion of concentration with a probability of 40% and 30% 
in the intermediate and high redshift bins, respectively. 
IConsehcc ( 2003.) argues that the concentration index re- 
veals the past star-formation history of a galaxy because 
of the correlation of bulge to total light ratios and stel- 
lar masses with this index. This behavior implies that 
the MlPS-non-detected galaxies have similar stellar pop- 
ulations and distributions as the LIRGs except for the 
LIRGs having a large population of massive young stars 
due to the presence of a starburst. 

4.2.2. Merging Galaxies 

It is also of interest to determine how many LIRGs are 
merging. We classify an object as a m erger if its asym - 
metry parameter is greater than 0.35 ijGonselicd l2003j) . 
The merger as defined here means a major merger, i.e., 
that the two progenitor galaxies have comparable masses. 
However, other strong asymmetries may also satisfy this 
numerical definition. By this definition, the LIRG sample 
is composed of 40±20 % merging systems, higher than 
the merging fraction (^15%) of the MlPS-non-detected 
sample. Although the merging fraction obtained by the 
asymmetry technique depends on the A value used for 
merger identification, this conclusion does not change 
with modest changes in this defining value of A. 

Our sample shows 40% of LIRGs at 0.3<z<1.4 are 
merging systems as identified by the CA method (which 
emphasizes high levels of morphological disturbance in 
systems with high surface brightness). This result 
is probab l y consi stent with but higher than that of 
iBell et al.l l|2005a|) who found that less than 30% of 
1500 star-forming galaxies at redshift 0.65<z<0.75 are 
strongly interacting based on visual classifications, and 
that of W. Zheng et al. (2004) who found that 6 (-17%) 
of 36 distant (z>0.4) LIR Gs in the CFRS are o bvi- 
ous ongomg major mergers. IW. Zheng et~an l|200l do 
not count the irregular galaxies (22%), which may have 
A>0.35, as merging galaxies. Another cause for these 
differences is that the morphologies of a merging galaxy 
at different surface brightness levels show different levels 
of asymmetry and thus visual identification on a rela- 
tively shallow image may mis-classify it as a non-merging 
galaxy. 

5. DISCUSSION 



5.1. Continuous Morphological Transformation plus 
Episodic Starbursts 

Our study of morphologies and infrared properties of 
high-redshift galaxies in the UDF has two basic results. 
The first one is that LIRGs and MlPS-non-detected 
galaxies in the UDF are predominantly late-type asym- 
metric galaxies. The LIRGs are somewhat more asym- 
metric on average than MlPS-non-detected galaxies and 
they are characterized by a higher fraction of very asym- 
metric {A > 0.35) systems and a lower fraction of sym- 
metric {A < 0.2) objects. This result is consistent with 
the proposal that non-symmetric structures play an im- 
portant role in enhancing the star formation and hence 
elevate the infrared luminosities of galaxies. The sec- 
ond result is that there is no one-to-one correlation be- 
tween the occurrence of high infrared luminosities and 
disturbed morphologies: highly asymmetric galaxies, in- 
cluding merging systems, do not necessarily have high 
infrared luminosities. 

On average, high-z LIRGs and normal galaxies are 
both more asymmetric than local irregular types, and 
their concentration indices are similar to the local late- 
type disk galaxies. Are the asymmetric normal galaxies 
and the LIRGs two different populations or are they the 
same population but in different stages? That is, are the 
activities of starbursts episodic so that the LIRGs are 
now in active stages and asymmetric non-LIRGs are simi- 
lar galaxies that are in quiescent states between-starburst 
stages? 

Based on the statistics at 0.3<z<0.9 listed in table 
2, 25% of asymmetric galaxies (A>0.2) are LIRGs at 
0.3<z<0.9. This fraction is ~10% in t he local uni- 
verse , assuming a merging fraction of 4% l|M arzke et alJ 
ITqQSI) . a LIRG fraction of 0.5% (Soifer c t al.lll986j) aid 
an asymmetric galaxy fraction in LIRGs of 74% based 
on the f raction of close pair and strongly interacting 
galaxies l|Sanders et al.llT988|l . One candidate popula- 
tion for asymmetric galaxies without starbursts is in- 
teracting galaxi es with low me tallicity and hence low 
dust. However, iKobulnickv fc KewlevI p004f) find that 
the shift of the luminosity-metallicity relation is smaller 
than 27% at z^l.O compared to the local trend, based 
on the nebular emission lines of galaxies (18.5 < Mb < 
-21.5) with SFR -0.1-10 Mp,yr-^. Th is result is consis- 
tent with previous stu dies l|Kobulnickv fc Zaritskvlll99^ 
ICarollo fc Lilhll200l . We conclude that both LIRGs 
and MlPS-non-detected galaxies are not likely to be 
extremely deficient of dust. Another candidate pop- 
ulation for asymmetric non-detected galaxies is inter- 
acting galaxies without gas, e.g., interactions between 
galaxies composed entirely of stars. However, such "dry 
mergers" appear in general to be uncommon 1 %) 
l)Bell et al.ll2005h() . Therefore, the similar distributions 
of concentration index for LIRGs and asymmetric non- 
MlPS-detected galaxies suggests that these two samples 
have similar formation histories. The moderately higher 
asymmetry for the LIRGs could arise from luminous off- 
nuclear star formation, or from extinction. Thus, the el- 
evated periods of star formation are probably occurring 
episodically. 

In the local universe, starbursts in a given galaxy 
have a short durat i on of 0. 01 to 0.07 Gvr llKrabbe et alJ 
IT99I iLutz et al.l IT99fit lAlonso-Herrero et alJ YimT . 
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Near-IR imaging and spectroscopy in the star- 
burst galaxies M82, IC 342 and NGC 253 indi- 
cate several episodes of enhanced st a r formation 
llRieke et alJ Il993t iSatyapal et all Il997t iBoker et all 
119971 lEngelbracht et alJ 11998': "^^ra tCT Schreib'CTn200n| r^ 
The anticorrelation between equivalent wi dth of Ha and 
galaxy pair spatial separation found by iBarton et al.l 
l)200Cl|) indicates that the starburst episodes are initiated 
by a series of close passes. In the early Universe, Ly- 
man Break Galaxies (LBGs) show comparabl e durations 
of starburst activi ty (Papovich et al. 2001 : Sha plev et al. 
1^05). Recently, lYan etraO (120041) and iPapovic h et al. 
l|2005i) found that the broadband photometry data for 
massive galaxies at Z'^1-3 are better fit by a model 
star-formation history with recent bursts superimposed 
on an underlying stellar population compared to sim- 
ple, monotonically-evolving stellar populations. They 
conclude that these galaxies have a complicated his- 
tory with many episodes of star formation. Theoreti- 
cally, simulations show that in merging galaxies, strong 
gas inflow is triggered and nuclear star formation is ig- 
nited only after enoug h gas is accumulated in the cen- 
ter of the g alaxies re .g.lNoguchilll991l:lMihos et alJll993l: 
iBarnes ~ Hcrnauisl ll99q) . During the whole merging 
process, strong enhancement of the star formation only 
occurs with a short duration, which supports the pro- 
posal for episodic starbursts. A star formation his- 
tory involving episodic events i s also indicated by re- 
cent semi-analytic st udies (e.g. iSomerville et al.l 1200 H 
iNagamine et al.ll2nn5j) . 

Given the merging timescale Emerge , the fraction / of 
Emerge at which galaxics are actually involved in star- 
bursts, and the timescale for episodic starbursts tESB, 
we can estimate the total number of episodes during one 
merging event N-ESB=f x tmorgc/^ESB- We adopt 0.5 Gyr 
for the timescale of the greatest merger activity (e.g. 
iPatton et a l."2000l . The characteristic timescale of star- 
bursts is estimated at 10-70 Myr as shown above. The 
fraction / can be computed from the merging fractions 
of LIRGs and infrared-inactive galaxies as listed in Table 
1 and the fraction of LIRGs in a given magnitude- limited 
sample. We adopt 15% as the fraction of LIRGs as shown 
in Section 4.1 and find / = 0.15 x 0.4/(0.15 x 0.4-^0.85 x 
0.16) — 0.3. Therefore, the average number of episodes 
during a merging event is around 4 for ^esb = 0.04 Gyr. 
This indicates that during each merging event, at least 
one episode of vigorous star formation can be triggered 
and multiple episodes of such activity are possible. 

This conclusion is supported by the studies of nearby 

galaxies cited above^ ^In addition, simulations by 

|Mihos & Hcrnquist 1)199(11) show that during merging 
events, bulgeless disc galaxies can trigger gas inflow even 
before the occurence of the actual collision of the galax- 
ies and thus two starbur st episodes can occur during the 
merging event. Recently. iTissera et al.l l^2002^ has shown 
that in hierarchial clustering scenarios, some merging 
events can have two episodes of starbursts: one owing to 
the gas inflow driven as the satellite approaches, and the 
second one occuring in the collisions of baryonic clumps. 

5.2. Comparisions to the Local LIRGs 

In the last column of Table 2, we present the statistical 
results f or local LIRGs list ed in the IRAS Bright Galaxy 
Sample ijSoifer et al.lll98^ . We extracted the image for 



each source from the Digitized Sky Survey (DSS) at the 
P0SS2/UKSTU blue band. The table shows that the 
morphologies of local LIRGs are similar to those of the 
MlPS-non-detected sample and the LIRGs in the UDF 
in terms of the fraction of late- type galaxies, of merging 
galaxies and of symmetric galaxies. The mean asymme- 
try and concentration of LIRGs in the UDF are slightly 
higher than those of local LIRGs. The K-S test shows 
that the probability that the LIRGs at 0.3<z<0.9 and 
local LIRGs have the same distribution of asymmetry 
and concentration are 3% and 9%, respectively. That 
is, the apparent differences are only marginally signifl- 
cant. In the local universe, a very low (4%) fraction of 
the total population is irregular and interacting galax- 
ies (jMarzke et al.l ll998^. As a result, the morphologies 
of local LIRGs are quite different from the local normal 
infrared-inactive galaxies, which is not the case at high- 
redshift. The small difference in the morphologies be- 
tween LIRGs and infrared-inactive galaxies in the UDF is 
due to the fact that, at high redshift, a significant fraction 
of galaxies are asymmetric systems in which starbursts 
are much more easily triggered. Locally, LIRGs are ex- 
ceptional, and an exceptional disturbance is required to 
trigger their elevated rates of star formation. Based on 
this trend, we can expect that at higher redshift (z>1.5), 
the morphologies of LIRGs and non-LIRGs will be even 
more indistinguishable. 

The similarity between the morphologies of local and 
high-z LIRGs suggests that a certain level of asymmetry 
is correlated with vigorous star formation, independent of 
redshift. This situation probably arises because once star 
formation is initiated, it evolves w ithin a galaxy inde - 
pendent of the galaxy environment l|Balogh et al.ll2?)03) . 
Therefore, a similar level of perturbation leads to simi- 
lar consequences in terms of elevated star formation. At 
z 1, most galaxies appear to have an appropriate level 
of asymmetry potentially to be transformed into LIRGs 
by relatively minor mergers. Therefore, the incidence of 
LIRGs is high, leading to a high density of star forma- 
tion. Locally, most galaxies are highly symmetric and 
only exceptional circumstances, such as an interaction 
with another massive galaxy, can produce the situations 
that are a prerequisite for a LIRG. Hence, the incidence 
of LIRGs is low and so is the density of star formation. 

6. CONCLUSIONS 

In this paper, we have investigated the morphologies of 
18 high-redshift (0.3< z <1.4) LIRGs in the UDF down 
to a rest-frame B-band surface brightness of hbo — 24.5. 
They are a representative sample of the LIRGs that dom- 
inate the star formation density at z ^ 1. We compare 
their morphologies to those of MlPS-non-detected galax- 
ies at similar redshift and of 49 local LIRGs. We used 
the concentration and asymmetry indices to quantify the 
morphologies of all these galaxies. Our main results are: 

(1) The fraction of LIRGs for Mb <-18.5 galaxies at 
this redshift range is around 15%, much higher than the 
0.5% in the local universe. Their morphologies are dom- 
inated by late-type galaxies with asymmetric (A>0.2) 
structures and include a significant number that appear 
to be mergers (40±20%). The morphologies of LIRGs 
in the UDF are similar to those of local LIRGs but pos- 
sibly somewhat more asymmetric. Optically luminous 
MlPS-non-detected galaxies at Z'^1 are nearly as asym- 
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metric as LIRGs, in contrast to the local universe where 
infrared-inactive galaxies are far more symmetric than 
local LIRGs. 

(2) All of the high-z galaxy types are relatively asym- 
metric, more so than even local irregular types. Their 
concentration indices are similiar to those of local late- 
type disk galaxies. We argue that the asymmetric LIRGs 
and asymmetric MlPS-non-detected galaxies are from 
the same parent population but in different stages of star 
formation, which suggests the star formation is episodic. 
Under this hypothesis, the fraction of the duration of 
the IR-active phase for asymmetric (A>0.2) galaxies is 
- 25%. 

(3) The morphologies of local LIRGs and LIRGs at 
z 1 are similar, suggesting that similar conditions 
within the galaxy lead to LIRG-levels of star forma- 
tion. At z ^ 1, such conditions arc common and LIRGs 



can develop by relatively minor mergers. Locally, most 
galaxies are highly symmetric, and only exceptional cir- 
cumstances such as an interaction with another massive 
galaxy can create the level of asymmetry associated with 
LIRG-levels of star formation. 
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Note. — Column (1): The MIPS ID of the LIRGs in the UDF. Column(2): RA of MIPS source. Column(3): DEC of MIPS source. Column(4): ACS ID of opitcal ^' 
counterpart of MIPS source. Column(5): The separation in arcsec between ACS source and MIPS source. Column(6): The redshift. Column(7): 24 flux. Column(8): e+ 
z-band magnitude. Column(9): Asymmetry. Column(lO): Concentration. Column(ll): The galaxy type of MIPS source based on CA system. Column 12: The photometric p 
redshift and the catalog number of X-ray counterparts of MIPS sources based on the CDFS 1 MS catalog by W. Zhcne ct al. |2003). Column 13: Where redshifts are from 
COMBO-17, the corresponding catalog number along with the magnitude at R band is listed; 'This work indi cates that the redshif t is computed based on our photo-z code 
iPere z-Gonzalez et al. 2005). 'Spec' means the spectroscopic re dshift from VLT/FORS2 spectroscopy catalog IVanzella et alj|2005ll and VIMOS VLT Deep Survey (VVDS) 
vl.O catalog liLc Fcvre et al...2004l) . 



TABLE 2 

The statistics of morphologies of Mb < —18.5 galaxies in the UDF and local LIRGs 
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Note. — ^Averages with standard errors of the mean. At the 0.3<z<0.9, both LIRG and MlPS-non-detected samples are 
complete. At the 0.9<z<1.4, both sample are incomplete. 
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Fig. 1.— The redshift distributions of LIRG (L/fl(8-1000/um)> IO^Lq), non-LIRG (Ljij (8-1000/im)< IO^^Lq) and MlPS-non-detected 
control samples. The dot-dashed histogram indicates the control sample, the dotted histogram indicates the LIRGs and the solid histogram 
indicates the non-LIRGs. The control and LIRG samples have almost the same redshift distributions, while non-LIRGs are mainly at low 
redshift due to the detection limit. The inserted plot shows the IR luminosity at the 24 nm detection limit (0.06 mjy ) as a function of 
the redshift. 
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Fig. 2. — The rest-frame B-band images (5x5 arcsec) of non-X-ray LIRGs in the UDF along with the redshifts, structural parameters 
and the galaxy types based on the CA classification scheme. The center of the image is the position of the MIPS source. The light solid 
and dashed curves show contours at 20% and 50% enclosing flux for IRAC and MIPS, respectively. The heavy solid and dashed circles 
show the position accuracy for IRAC sources (0.1 arcsec) and MIPS sources (0.5 arcsec). 
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Fig. 3. — The distributions of LIRGs (triangles), non-LIRGs (open eircles) and MlPS-non-detected control sample galaxies (stars) in 
the plane of concentration and asymmetry. The two solid curves divide the plane into three regions populated by early-, intermediate- and 
late-type galaxies. The dotted line indicates A = 0.35, to the right of which galaxies are identified as merging systems. The cross shows 
typical errors of concentration (0.3) and asymmetry (0.04). 



